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ABSTRACT 



Aims. We use 9.5-yr of BiSON Sun-as-a-star data to search for dependence of solar-cycle parameter changes on the angular degree, 
of the data. The nature of the Sun-as-a-star observations is such that for changes measured at fixed frequency, or for changes averaged 
across the same range in frequency, any { dependence present carries information on the latitudinal distribution of the agent (i.e., the 
activity) responsible for those changes. 

Methods. We split the 9.5-yr timeseries into contiguous 108-d pieces, and determine mean changes in the damping of, power in, and 
energy supplied to the modes through the solar cycle. We also apply a careful correction to account for the deleterious effects of the 
ground-based BiSON window function on the results. 

Results. From our full analysis we obtain a marginally significant result for the damping parameter, where the mean change is found 
to be weakest at ^ = 0. The other parameters show hints of some dependence in L 

Conclusions. Our main conclusion is that the mean fractional solar-cycle change in the ( = damping rates is approximately 50 % 
smaller than was previously assumed. It had been common practice to use an average over all low-f modes; our downward revision 
of the radial-mode value has implications for comparisons with models of the global solar cycle changes, which are usually based on 
a spherically symmetric geometry. 

Key words. Methods: data analysis - Sun: helioseismology - Sun: activity 



1. Introduction 

The fact that damping rates and powers of the global p 



(e.g., iQiaolinetal., 2000 


Kommetal., 2000; Salabert et all 


20031 Ijimenez-Reves et al. 


2003; Jimenez-Reves et al..> ,20041 


Salabert & Jimenez-Revesl 


2006|). 


At least where the main part 



of the p-mode spectrum is concerned damping gets heavier, 
and observed power gets weaker, as the level of solar activ- 
ity increases. Information on the damping and power parame- 
ters comes straightforwardly from the observations. The damp- 
ing rates are assumed to be linearly related to the linewidths of 
the resonant peaks in the frequency power spectrum; while the 
powers are proportional to the product of the peak widths and 
heights (the latter more formally termed the maximum power 
spectral densities)^ 

[Komm et aT] (l2002t ) took advantage of the large number of 
components available to analyze in the medium-degree range 
(from { = 40 to 80) the latitudinal dependence of the width and 
height changes. They found that the changes were concentrated 
in latitudes occupied by the active regions. The damping rates 
and powers, like the mode frequencies, therefore seemed to be 
responding in some fashion to changes wrought on the active 
regions by the changing magnetic fields. 

Most of the low-^ results have come from the Sun-as-a-star 
data (e.g., from BiSON, GOLF, IRIS and VIRGO/SPM; or from 
GONG and MDI data combined into a Sun-as-a-star-like proxy). 
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These data show only the even t+m components; the odd compo- 
nents are so weak as to be unobservable. Since it is the outer, sec- 
toral components (with { = \m\) that appear most prominently, 
the Sun-as-a-star mode parameters estimated by the usual anal- 
ysis methods are dominated by, and therefore close to, the sec- 
toral values. This characteristic gives a noticeable change in the 
latitudinal sensitivity from { = Oto^ = 1,2 or 3. The spatial 
sensitivity of Sun-as-a-star data at the latter three values of { is 
weighted toward the lower latitudes, where the active regions re- 
side. The Sun-as-a-star parameters at these { are therefore more 
sensitive to the solar cycle than are the £ - data. Evidence 
for t dependence in the shifts of the Sun-as-a-star mode param- 
eters therefore carries information on the latitudinal distribution 
of the agent responsible for those shifts. (Changes in inertia, at 
fixed frequency, between these £ are very small indeed.) 

It is much more difficult to uncover spatial dependence of the 
mode parameter changes in the low- than in the medium-^ data 
because there are far fewer components to analyze, and uncer- 
tainties on the results are commensurately larger When changes 
to damping and power were first uncovered in the Sun-as-a-star 
data, results were averaged over { to reduce errors. Given the 
modest precision in the results, it was assumed the ^-averaged 
values also provided a working proxy of the radial-mode shifts. 
And so these averages were used as meaningful c omparisons for 
models of global changes in damping (see, e.g.. lHoudeket all 
2001), models set up in spherical geometry, i.e., pertinent only 
to the { - case. With better analysis, and more data, it has 
become possible to give results for each {, and to therefore test 
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whether the £ - shifts really are weaker than an average across 
^ = to 2 or 3. 

We made a first cut at such an analysis in IChaphn et alJ 
(120031 . Our results, from BiSON Sun-as-a-star data, suggested 
very strongly that the f = linewidth changes were indeed sig- 
nificantly weaker than at{ - 1 and 2. The clear implication was 
that results of the theoretical models would now need to be com- 
pared to this new, smaller shift. In the analysis of the BiSON 
data, we had to allow for the corrupting influence on the results 
of the ground-based BiSON window function. We did so by ap- 
plying corrections that were designed originally for a different 
study, using data from a diff'erent period. In this paper we revisit 
our analysis. We design and implement a correction procedure 
for the 9.5-yr BiSON dataset in question. We find that implemen- 
tation of this internally consistent, and more accurate, correction 
gives little change to the results, and reinforces our earlier con- 
clusion. In summary, the mean fractional solar-cycle change in 
the £ - damping rates is approximately 50 % smaller than was 
previously assumed. 

2. Data analysis 

2.1. Observations 

Our analysis made use of 3456 d of data on the low-angular- 
degree solar p-modes, collected by the six-station, ground-based 
Birmingham Solar-Oscillations Network (BiSON). This 9.5-yr 
dataset had a starting date of 1992 July 26, and spanned the 
falling phase of solar cycle 22 and the rising phase, and the max- 
imum, of solar cycle 23. 

The instruments at each BiSON site make Sun-as-a-star ob- 
servations of the Doppler shift of the potassium Fraunhofer line 
at 770 nm (e.g., .Chaplin et al.. 1996). Ra w data were fi rst pro- 
cessed in the manner described by Elswort h et alJ ( Il995h to yield 
daily calibrated velocity residuals. We then combined coherently 
the resulting a; 2 x 10"* individual daily sets from all six stations 
to yield the principal 3456-d time series of residuals. The duty 
cycle of this combined set was ~ 76 %, with breaks in coverage 
largely the result of inclement weather 

In order to derive solar-cycle trends, the 3456-d series was 
divided in thirty-two contiguous, independent 108-day pieces, 
which had a range of duty cycles between 63 and 87 per cent. 
This 108-d interval, a multiple of the solar rotation period, was 
chosen to minimize effects dependent on the solar rotation. 



2.2. Mode parameter extraction 

The power spectrum of each 108-day time series was fitted to 
yield estimates of the mode parameters over a range in frequency 
1600 < V < 3800 yuHz. This fitting was do ne by a multi-step iter- 
ative method base d on lGelly et aP (^2002^ . The asymmetric pro- 
file of E^iaam & Kosovich'enT998l) was used to describe each 
component, as: 



Mix) = H 



(1 +bx)^+b^ 

1 -Hx2 ' 



(1) 



where x = 2(v - vo)/r, and H, F, vq are the mode height (max- 
imum power spectral density), the fwhm linewidth, and the cen- 
tral frequency respectively. The parameter b is the asymmetry 
coefficient, which in the present analysis was fixed at small neg- 
ative values determined by analysis of longer datasets. 

Due to their close proximity in frequency, modes were fit- 
ted in pairs (i.e., £ - 2 with and ( - 1> with 1). An addi- 
tional offset was included in the fitting model to describe the 



background level in the fitting window. The first temporal side- 
bands at 1 1 .57 jt/Hz were also included in the model. The mode 
parameters were extracted by maximizing a likelihood function 
commensurate with the 2 degrees-of-freedom statistics of 
the power spectrum. The natural logarithm of the height, width 
and background terms were varied - not the parameter values 
themselves - in order to give quasi-normal fitting distributions. 
Formal uncertainties on the fitted values were then derived from 
the Hessian matrix of each fit in the usual manner. 



2.3. Damping and excitation parameters 

A generally accepted analogy for the p modes is a forced and 
damped harmonic oscillator For any given mode in, {), the mode 
width, r„/, is then linearly related to its damping rate, ri„ (. The 
total velocity power, P„/, gives a measure of the balance be- 
tween the excitation and damping, and is proportional to the 
product of the height, //„ f, and the mode width, F,, The rate at 
which the energy is supplied to the modes, E„j, is a direct mea- 
sure of the net forcing and is independent of the mode damp- 
ing. This 'supply rate' can be computed as the product of the 
mode width, F,,/, and the mode energy, E„^e, the mode energy 
being proportional to the mode power, f „ For a detailed de- 
scrip tion of the mo de damping and excitation parameters, see, 
e.g., Chaplin et al. (2000). 

Since the peak-finding procedure returned natural logarithms 
of the mode heights and widths, estimates of the mode velocity 
power and supply rate were made by taking the combinations: 



log,(P„/) = log,(//„,f) + log,(F„,f) -H Ci, 
and 

logg(£'„/) = log^X^^«,f) + 21og^(r„,^) + C2. 



(2) 



(3) 



Variations in the fitted logarithmic parameters therefore cor- 
responded to fractional variations of the absolute parameters. 
The factors Ci and C2 in the above contain information on the 
mode inertia, and mode visibility given by the instrument. The 
factors were assumed to remain unchanged over the epoch in 
question, and were therefore ignored in the solar-cycle analy- 
sis. Uncertainties on the mode velocity power and on the mode 
energy supply rate were determined by taking into account the 
strong anti-correlation between the extracted mode heights and 
mode widths ( ChaoUn et al...2000) . 

Hereafter, we write the natural logarithm of H„ [, F„_^, P„/, 
and E„x as h„/, p„x, and e„x respectively. 

2.4. Determination of the fractional variations 

Estimates of the fractional variations of the mode parameters 
over time were obtained by performing a linear regression of 
each of the fitted parameters y„/, p„ (, and e„/, on the 10.7-cm 
solar radio flux' (Fioj). The radio flux was used as a robust 
proxy of the global level of solar activity. The best-fitting gra- 
dients therefore gave for each mode a measure of sensitivity of 
its parameters to the solar activity. 

The regressions were computed as both un-weighted and 
weighted linear fits. The formal parameter uncertainties from the 
mode-fitting procedures, derived from the Hessian matrices of 
the fits, were used here to compute weights for weighted linear 
regressions. Fitted gradients from the regressions were then aver- 
aged at each degree over the frequency range between 2600 yuHz 
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and 3600 /zHz. This averaging yielded estimates of mean mode 
changes along the cycle. 

There are different ways to compute these averages and their 
associated uncertainties. On the one hand, a simple, unweighted 
mean of the input residuals may be calculated, with the scat- 
ter on the residuals used to determine the uncertainty on the 
mean (in the usual manner). An alternative is to use the Hessian- 
determined uncertainties on the best-fitting gradients to calculate 
weighted means. The error on the mean can then be estimated 
in one of two ways: first, from a suitable combination of the 
Hessian-determined uncertainties, giving an internal error; and 
second, by using the Hessian-determined uncertainties to weight 
a calculation of the dispersion of the data, giving an external 
error (Chaplin et al., 1998). 

Since the regression fits were performed in two ways (un- 
weighted and weighted) there are six ways in all to estimate 
the degree dependence of the p-mode damping and excita- 
tion parameters, and the associated uncertainties (although the 
weighted-average methods of course give the same means, but 
have different errors). 



3. Correction for bias 

Gaps present in the timeseries of a ground-based window func- 
tion lead to overestimation of the mode widths and underesti- 
mation of the mode heights dChaplin et al .. 2003b). Proper al- 
lowance must therefore be made for any bias introduced by the 
window functions of the thirty-two 108-d BiSON timeseries. 
Furthermore, marked changes in the intrinsic quality of the data 
can in principle influence the parameters, via the impact changes 
in quality may have on the observed height-to-background ratios 
of the mode peaks. Changes in height-to-background ratio must, 
however, be large to give any significant bias in the parameters. 

We have used Monte Carlo simulations of artificial Sun-as-a- 
star data to help us correct the BiSON results for bias introduced 
by these effects. We found that the bias was dominated by the 
impact of the different 108-d window functions; the impact of 
changes in data quality was far less significant. The domination 
of the window function is illustrated in Fig.^ It shows the mean 
background levels, of the thirty-two 108-d BiSON, spectra over 
the range 2600 to 3600 /vHz (-■-). The values have been scaled 
so that the highest mean background level is set to unity. The 
levels are plotted as a function of the duty cycle, D, of each of the 
sets. As fill decreases, the redistribution of power by the window 
function increases the background level. 

Also shown in Fig. ^ are results for a single, 108-d GOLF 
timeseries. The GOLF timeseries has nearly a 100% duty cy- 
cle, and has been modulated, in turn, by each of the thirty-two 
BiSON window functions. The resulting mean background lev- 
els have been plotted (-o-). Because the GOLF data are the same 
for each computation, the change in the background level comes 
wholly from the impact of the window function. Changes to the 
background in these modulated GOLF data are seen to follow 
very closely the changes observed in the BiSON data. From 
this, we conclude the window function dominates changes to 
the background. The result simplifies the generation of artificial 
data for the Monte Carlo simulations in that we do not need to 
worry about introducing changes in intrinsic data quality across 
the simulated datasets. 



BiSON 

Sr- GOLF with BiSON window functions 




0.6 0.65 0.7 0.75 0.8 0.85 0.9 

Fractional fill D 

Fig. 1. Normalized background level for the thirty-two 108-d 
BiSON spectra (-■-) and the single 108-d GOLF set modu- 
lated in turn by each of the BiSON window functions (-o-). 
Background levels are averaged over 2600 and 3600 juHz, and 
normalized so that the highest datum takes the value unity. Data 
are plotted as a function of the duty cycle, D. 



3.1. Bias estimation from artificial Sun-as-a-star data 

The artificial data were generated using the solarFLAG^ 
(Chaplin et al., 2006) mode simulation code. The Laplace trans- 
form solution of the equation of a forced, damped harmonic os- 
cillator was used to generate each component at a 40-s cadence 
in the time domain, in the manner described by IChaplin et alJ 
( 1997). Components were re-excited independently at each sam- 
ple with small 'kicks' drawn from a Gaussian distribution. The 
model gave rise to peaks in the frequency power spectrum whose 
underlying shapes were Lorentzian. It should be pointed out that 
the observed low-^ peaks are slightly asymmetric in shape (al- 
beit at the level of only a few per cent at most). 

Datasets were constructed component by component in the 
time domain. Sets were made with a full cohort of simulated 
low-^ modes, covering the ranges < f < 5 and 1000 < v < 
5500 juHz. A database of mode frequency, power and linewidth 
estimates, obtained from analyses of GOLF and BiSON data, 
was used to guide the choice of input values for time-series con- 
struction. A rotational frequency splitting of 0.4yuHz was im- 
posed on all the non-radial modes, irrespective of overtone num- 
ber and angular degree, to match that extracted from observa- 
tions. A 'pink' background noise component was also added in 
the time domain, the size of which increased at lower frequen- 
cies in order to give signal-to-background ratios commensurate 
with GOLF or BiSON-like Doppler velocity data. 

Two cohorts of data were made, both comprising twenty 
3456-d timeseries. Each timeseries was excited with indepen- 
dent realization noise. We introduced into the first cohort (#1) 
systematic changes in the frequency, damping and power pa- 
rameters so as to si mulate the effects of the solar cycle (see 
IChaplin et all l2006h . These changes were based on variations 
of a modelled proxy of the 10.7-cm radio flux. The proxy was 
made to mimic the radio flux for the real epoch covered by 
the BiSON timeseries. Parameter variations were calibrated to 
give equal, and oppositely directed, changes in the damping and 

^ FLAG URL: http ://bison.ph.bham.ac.uk/i ~ wjc/Research/FLAG.html 
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power param eters (i.e., 6y„, { = 6p „,(), to match previous obser- 
vations (e.g.. IChaplin et all l2000h . In this scenario the supply 
rate, e„ /, shows no net change. The simulated changes were fre- 
quency dependent, showing their largest variations at 2900 yuHz. 
Furthermore, they were given no t dependence. In contrast, in 
the second 'stationary' cohort (#2) all mode parameters were 
held constant in simulated time. 

We analyzed each artificial 3456-d dataset in the manner de- 
scribed in Sections [2.21 and the first step being to split them 
into 108-d pieces. From the analysis, we obtained ^-dependent 
estimates of mean changes in the y„/, p„ (, and e„/ parameters. 
Results from all 3456-d sets were then averaged within their re- 
spective cohort to reduce errors. Final estimates of the mean pa- 
rameter changes could then be compared with the true changes 
introduced, to allow for a determination of the bias. 

Similar mean values and errors were returned by the differ- 
ent averaging methods (Section [2. 4> . For the sake of clarity we 
show in what follows results from the weighted regressions only. 
The three panels in the left-hand column of Fig.|3show the mean 
estimated parameter variations, per unit change in the radio flux, 
given by analysis of artificial cohort# 1 (the 'solar-cycle' data). 
The right-hand panels show results from cohort #2 (the 'station- 
ary' data). The horizontal dashed lines show the actual mean 
changes that were introduced in the artificial parameters. The 
plotted uncertainties in each panel show the error bars on the 
mean results of the twenty artificial datasets. 

The results in Fig.|3show evidence for both overestimation 
and underestimation of input values, dependent on the parame- 
ter and cohort. Bias is most severe in the determination of the 
change in power and supply rate in cohort#l. In spite of the 
off'sets, it is encouraging to note that there is no appreciable vari- 
ation of the bias with degree. Any variation would certainly be 
well within the uncertainties expected from analysis of a single, 
3456-d dataset (which are a factor V20 larger than the errors 
plotted in Fig.|5J. 

From these artificial data we then went on to derive a 
window-function parameter correction that could be applied to 
the real BiSON observations. Because of the small errors on the 
artificial data results, the correction would be well constrained, 
and uncertainties on the results in Fig.|2would therefore not add 
significantly to the uncertainties on the final, corrected BiSON 
results. 



3.2. Formulation of correction procedure 

The following procedure was applied to the fitted parameters 
within each cohort of artificial data. 

We began by computing the ratio of the fitted and input pa- 
rameters of every fitted mode. Ratios were made for the heights, 
h, and widths, y. Recall that each of the twenty independent 
3456-d sets was split into thirty-two 108-d pieces for fitting. 
We therefore obtained many ratios, Qh(i, j, n, () and Qy{i, j, n, {). 
The variable / ran over the independent sets (1 < / < 20); while 
j ran over the 108-d pieces (1 < /' < 32), each having their own 
distinct BiSON window function. 

Next, we averaged the ratios over the twenty independent 
datasets (over ;)■ This gave a total of thirty-two mean parame- 
ter ratios for every mode, one for each 108-d window function. 
These mean parameter ratios, Q{j, n, (), were then regressed lin- 
early against the duty cycles, D(i), of each of the windows; in 
fact, we fitted to 1 - D{i), rather than D{i). The fits gave best- 
fitting gradients 6Q^{n, {). 



By following the above procedure for each cohort of artificial 
data, we obtained two sets of mean ratios and best-fitting gradi- 
ents. We were then in a position to construct, from the ratios and 
gradients, a formalism for correcting the fitted BiSON parame- 
ters for the window function. In short, it would allow us to make 
estimates of the unbiased parameters by linear extrapolation to 
D = 1.0. 

First, we checked the procedure could be applied to the fit- 
ted artificial mode parameters to recover correct estimates of the 
mean input changes in cohorts #1 and #2. This would involve 
performing a full analysis on each of the artificial 3456-d sets, 
and would be a dry run of the full correction and analysis as it 
would be applied to real BiSON data. Taking the width param- 
eter as an example, estimates of the unbiased widths, jQij, n, €), 
were given by the linear combination: 

jQij, n, {) = y{ j, n, €) X {e,(j, n, ^) + [1 - D(/)] X c5e°(n, ^))-'.(4) 

A similar equation was used to define a correction for the mode 
heights. Once we had obtained estimates of the unbiased widths 
(and heights) from Equation|2(and its mode-height counterpart), 
the corrected parameters were then used to estimate mean pa- 
rameter variations. We used the corrected parameters as input to 
the analysis outlined in Section lT^ and found we were able to 
recover, within errors, the input cycle shifts when we analyzed 
datasets from cohort#l; and the expected null changes when 
we analyzed datasets from cohort#2. To illustrate the success 
of the correction procedure, Table^shows the linear correlation 
between the estimated and input width changes for sets in co- 
hort #1 before (top row) and after (bottom row) the correction 
was used. Given the numbers of data involved, these represent 
significant improvements in the correlation, and the final result. 

Table 1. Linear correlations between the fractional variations in 
mode width (% per unit of change in F10.7) extracted from the ar- 
tificial datasets in cohort# 1 and the true, input variations without 
and with temporal window function corrections (Equation|3. 





t = Q 


I = 1 


e = 2 


( = 3 


No correction 


0.80 


0.93 


0.96 


0.61 


With correction 


0.98 


0.99 


0.99 


0.94 



4. Results with Sun-as-a-star BiSON data 

With the correction checked on the artificial data, we proceeded 
to apply the procedure to fitted parameters from the thirty-two 
108-d BiSON datasets. We used the mean parameter ratios, 
Qij, n, (), and the best-fitting gradients, 6Q°{n, £), from the ar- 
tificial data. Recall there were two sets of these coefficients: 
one from artificial cohort# 1, the other from artificial cohort#2. 
The window-function correction procedure was therefore ap- 
plied twice, yielding two sets of corrected parameters. 

Estimates of the unbiased height and width parameters were 
then analyzed, as per Section 12.41 This give the sought-for es- 
timates of the mean parameter changes, per unit change in the 
radio flux. Again, these changes - which are shown in Fig.|3- 
represent averages for each ( over the range 2600 to 3600 //Hz. 
The plots in the left-hand column show results given by the co- 
hort #1 correction coefficients; those in the right-hand panels 
results given by the cohort#2 coefficients. The results are also 
summarized in Table |2] Here, we have multiplied the shift-per- 
unit-activity values shown in Fig.|3(the external mean data -♦-) 
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by the change in the radio flux observed across the thirty-two 
108-d pieces, to give estimates of the mean parameter changes 
over the solar cycle. 

We found good agreement, to within the errors, between 
results given by the cohort #1 and #2 corrections. Averages 
made over the I - 0, 1 and 2 results were in reasonable 
agreement with previous analyses of the BiSON, IRIS and 
GOLF Sun-as-a-star data (Chaclin et al., 2000; Salabert et al., 
[2003; Jimenez-Reyes et al., 2003; Jimenez-Reyes et al., 2004). 
However, the most striking feature of Fig. |3]is the apparent rel- 
ative weakness of changes in the { = Q widths, implying an { 
dependence of the results. Inspection of results for the power pa- 
rameter showed hints of a mirrored, albeit weaker, trend. There 
were also indications of a jump in { in the supply-rate results (in 
particular for the cohort # 1 corrected data). 

To ascribe some formal significance to the weakness of the 
£ - Q width shift, we first compared its value with the weighted 
mean of the f - 1 and 2 width shifts (we ignored the less-well- 
constrained { = 3 data). The { - value was found to be about 
2a- different from the mean over f = 1 and 2, for both cohort cor- 
rections. However, when we computed a linear regression of the 
mean shifts on £, we found gradients x 3cr and x Act different 
from zero for the cohort #1 and cohort #2 coefficients respec- 
tively. We conclude that there is marginally significant evidence 
for an € dependence of the width changes. Similar analysis of the 
other parameters gave weaker levels of correlation, and therefore 
showed no evidence for significant variation, with £. Whether 
the hint of variation in the power and supply rate results is real 
only more data, and independent confirmation from other obser- 
vations (e.g., GOLF) will tell. 

5. Conclusion and discussion 

We have analyzed some 9.5-yr of BiSON Sun-as-a-star data - 
collected in solar cycles 22 and 23 - to search for dependence 
of the mode excitation and damping parameter changes on the 
angular degree, £, of the data. The nature of the Sun-as-a-star ob- 
servations is such that for changes measured at fixed frequency, 
or for changes averaged across the same range in frequency, any 
£ dependence present carries information on the latitudinal dis- 
tribution of the agent (i.e., the activity) responsible for those 
changes. 

We split the 9.5-yr timeseries into contiguous 108-d pieces, 
and determined mean changes in the damping of, power in, and 
energy supplied to, the modes through the solar cycle. We also 
applied a careful correction to account for the deleterious effects 
of the ground-based BiSON window function on the results. This 
correction was calibrated by, and then fully tested on, artificial 
seismic data generated by the solarFLAG mode simulation code. 

From our full analysis we obtained a marginally significant 
result for the damping parameter, where the mean change was 
found to be weakest at f = 0, and higher in data on the other 
£. The result implies the damping is strongest in the active re- 
gions, confirming r esults on the more numerous higher-/" data 
iKomm et al.Ll2002l) . The other excitation and damping parame- 
ters we investigated showed hints of some dependence in £, but 
nothing that could judged as statistically significant. 

Our main conclusion is that the mean fractional solar- 
cycle change in the £ - Q damping rates is approxi- 
mately 50 % smaller than was previously assumed. It had been 
common practice to use an average over all low-/ modes 
(where mean solar-cycle shift values have averaged abo ut 18% 
(ChapHn et al., 2000; Salabert et al., 2003; Jimenez-Reve s et all 
|2003>: Jimenez-Reves et al.> .2004) . Our downward revision of 



the radial-mode value has implications for comparisons with 
models of the global solar cycle changes, which are usually 
based on a spherically symmetric geometry. 
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Table 2. Fractional variations (%) of the low-degree p-mode damping and excitation parameters from solar minimum to solar 
maximum for each I value, averaged over 2600 and 3600 //Hz. Final results are shown for application of the window-function 
correction coefficients from artificial data cohorts # 1 and #2 (see Sec. l3.2t . Data correspond to the external mean shifts, rendered 
as -♦- symbols in Fig.|3 



Parameter 1 = t=\ t = 2 ^ = 3 <f = 0, 1,2) 



Correction coefficients from cohort # 1 

5y 9.9 ±3.0% 22.5 ±5.2% 25.4 ±6.7% 20.9 ±13.3% 14.7 ± 2.4% 

5p -8.7 ±4.7% -10.9 ±3.7% -16.0 ±3.2% -10.7 ±4.8% -12.7 ±2.2% 

6e -1.3 ±3.8% 10.6 ±4.3% 9.5 ±5.1% 9.6 ±12.7% 5.2 ±2.5% 

Correction coefficients from cohort # 2 



8y ^.\±'iA% 18.6 ±5.5% 22.6 ± 6.0% 16.4 ±13.2% 13.8 ±2.6% 

8p -^A±'i.'i% -1.1±2,.6% -n.4±3J% -13.2 ±4.7% -11.3 ±2.0% 
6e -4.9 ±3.1% 8.8 ±4.8% 4.8 ±6.4% 1.7 ±14.2% -OA ±2.4% 
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Fig. 2. Results of analysis of artificial data. Shown are mean fractional variations (% per unit of change in F10.7) in mode damping, 
6y (upper panels), mode velocity power, dp (middle panels), and mode energy supply rate, 6e (lower panels), as a function of the 
mode angular degree {. Panels in the left-hand column show results averaged over the twenty artificial 3456-d sets in cohort # 1; 
those in the right-hand column similar averages over sets in cohort#2. The dashed lines in each plot represent the modelled, input 
changes. The symbols give information on how the averages were computed (see Sec. l2.4> : unweighted mean -o-; weighted mean 
with internal error weighted mean with external error Plotted errors are those on the mean of the twenty artificial datasets. 
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Fig. 3. Results of analysis of BiSON data. Shown are mean fractional variations (% per unit of change in Fioj) in mode damping, 
6y (upper panels), mode velocity power, dp (middle panels), and mode energy supply rate, 6e (lower panels), as a function of the 
mode angular degree {. Panels in the left-hand column show results given using the cohort # 1 coefficients in the window-function 
correction; those in the right-hand column results from using the coefficients from cohort#2. The solid horizontal lines in each 
plot mark the zero-change level for reference. The symbols give information on how the averages were computed (see Sec. l2.4> : 
unweighted mean -o-; weighted mean with internal error weighted mean with external error 



